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0^ , A theory accounting for the specific features of magnetic polarons (MP) in the presence of spin glass 

^~^ • order is presented. We derive and solve self consistent equations for i) the polaron magnetisation, 

C^ ' ii) the thermodynamically averaged carrier-spin, and iii) for the spin glass order parameter. The 

D , temperature dependence of these quantities is analysed in detail. The modification of the spin 

pL^ ' glass phase due to the presence of the exchange field of the carrier inside the magnetic polaron 

-v^ I volume is investigated. The onset of spin glass order leads to a plateau-like flattening in the 

temperature dependence of the MP energy at low temperatures. It is found that solutions of spin 

I I ' glass equations are needed to optimally flt the experimental data of the temperature dependence of 

C , the exciton magnetic polaron (EMP) energy in (Cd,Mn)Te. Moreover, the dynamical aspects of the 

K . MP formation are discussed. Our model predicts qualitatively different temperature dependences 

tyj ' of the MP formation time in different dynamical scenarios. 
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Diluted magnetic semiconductors (DMS) exhibit a variety of magnetooptical effects (giant Zeeman splitting, giant 
' ^ Faraday rotation etc.) which allow for the investigation of these materials with both standard magnetometric and 
optical jteftiuiiques (for detailed review see Ref.El). The DMS (Cd,Mn)Te has been studied intensively during the last 
decadea3~ll3, its magnetic properties are determined by the spin 5/2 of the Mn^+-ions substituting Cd-ions randomly 
on the geometrically frustrated fee lattice. Excitons are prominent quasiparticles in this material and probe locally 
the magnetic state of the system. Under suitable conditions they align the surroundig manganese spins and turn into 
polarons. The properties of these new excitations are determined by a variety of sample parameters like quality and 
J^ , structure of barrier surfaces. However, in all probes with sufficiently high manganese content (or at temperatures 
,__! ' below 1 K) the disorder induced spin glass (SG) phase must be taken into account. We focus on this aspect of 
f^ , polaronic DMS physics. 

OQ ■ The magnetically active manganese ions are coupled via an antiferromagnetic interaction whose strength falls off 
^P ' rapidly with distance (the next-nearest neighbor (NNN) coupling is one order smaller than the nearest-neighbor 
(NN) interaction, which equals —5.7KB). By increasing the manganese concentration, clusters of magnetically coupled 
Mn-ions are formed. At concentrations above 20% a percolating Mn-cluster emerges, exhibiting a transition into a 
"t^ , spin glass phase. |a|l_low (helium) temperatures, which is caused by the combination of geometric frustration and 

a site randomnessLfSll3. The SG phase is characterised by the order parameter q = limt-,oo{S{0)S{t)), where t is the 

observation time, 5(0) and S{t) denote spins at identical sites and at times and t, and brackets represent the 
^-^ • thermodynaHiic average. In zero field q is zero in the paramagnetic phase and assumes a nonzero value in the spin 

Ch ', glass phasellj. 

O ' Experimentally the SG phase transition has been observed by a cusp in the temperature dependence of the magnetic 
susceptibility and by the onset of irreizEraihility in the field cooled and zero field cooled measurements of the remanent 
magnetisation (for reviews see Refs.ErEflo). The explicit use of spin glass equations for fitting experimental data of 
the magnetic polaron energy is however new and a key point of the present work. 

, . , The manganese ions interact strongly with electronic states in the valence and conduction bands. In magnetic fields 

j^ ' this leads to the giant Zeeman splitting of the band states, which extends up to tenths of meV. In the absence of an 
external magnetic field the formation of the magnetic polaron (MP) state - a carrier with a cloud of magnetically 
polarised Mn spins which has a lowered energy - becomes possible. The dominant contribution to this energy gain 
comes from the exchange interaction of manganese spins with the hole-spin. Magaetic polarons were investigated 
in detail for the case of Cjarriers bound to a shallow impurity (for review sec Rcf.l3) and for excitons localised by 
compositional fluctuationaiil. 

Since the properties of the magnetic polaron are strongly influenced by the magnetic state of the system, one can 
expect that the formation of the spin glass phase should affect the properties of this quasiparticle. The onset of SG 
order should show up for example in the temperature dependence of the MP energy. The aim of the present work 
is to analyse theoretically the effects related to the interaction between the magnetic polaron and SG order, proving 
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that SG effects are specific, hence distinguishable from any kind of antiferromagnetic clustering effects, and indeed 
required to explain qualitatively and quantitatively the new experimental data. 

The paper is organized as follows. In section || a model which accounts for the interaction between spin glass and 
magnetic polaron is presented and selfconsistent equations are derived for this model. An analytical analysis of these 



equations is presented in section III followed by detailed numerical evaluations and a comparison of the theoretical 
results with the experimental data in section |^ Finally in section ^ an extended model, which takes into account 
dynamical aspects of magnetic polaron formation in the presence of spin glass order is presented. 

II. THE MODEL AND ITS COUPLED POLARON SPIN GLASS EQUATIONS 

In this paper we-epasider the magnetic part of the MP energy, which is directly influenced by the SG ordering of the 
Mn spin systemOo. The exchange energy can also lead to the shrinking of the wave function of the polaron, thus 
contributing indirectly to the change of the kinetic energy of the MP. This contribution will be neglected tlpp|iia;tL0ut 
this paper. Furthermore, the contribution of the hole dominates the exchange energy of MP in (Cd,Mn)Teli3'E3. For 
simplicity we account for the hole contribution only, neglecting the small part of the exchange energy gained via 
coupling between electron and spin glass. These assumptions are surely justified for localised excitons. 
We adopt here the Villain model of Ising pseudospinaO'EI. In Ref.H detailed arguments were given in favor of this 
model and its applicability to (Cd,Mn)Te. The lower critical dimension for its Ising SG transition being smaller than 
three together with the fact that the Dzyaloshinskii-Moriya interaction, which would also turn the Heisenberg into 
an Ising SG (provided this coupling is isotropic itself; anisotropy would lead to special universality classesEl), is too 
small, supports the above model description. In the model considered here each pseudospin belongs to a tetrahedron 
with Mn spins on the corners. Due to the nearest neighbor antiferromagnetic interaction between the Mn spins the 
tetrahedron is frustrated and its ground states can be characterised by a mixed product ^ijk = Ui ■ (ctj x CTk) of three 
neighboring spins which can assume only two opposite values. The normalised value of this product S = ^/|^| can be 
1 or —1 and is hence considered as an Ising pseudospin. 

The interaction between hole- and Mn-spins will be described by the following model. Consider a percolating cluster 
with N pseudospins in (Cd,Mn)Te with random frustrated interaction between the spins leading to spin glass order 
at low temperatures. We assume the interaction between Ising pseudospins to be of infinite range. The hole spin of 
the exciton interacts with a part of iVi = aN pseudospins of the entire spin cluster. In this way the finite range of 
this exchange coupling is taken into account. We represent this Kondo-like coupling by —{k/Ni)s^ J2i2i Si- Here k, 
s^ , and Si denote the effective exchange constant, z-component of the hole spin, and the Ising pseudospin at site i 
respectively. The "exchange box" -model is assumed, i.e. the wave function of the hole is constant inside the magnetic 
polaron and zero outsideEll. Our assumption implies |^(r)p = 1/A^i inside the polaron and zero outside. Therefore 
the Hamiltonian of the model is given by 

ATi N 

if = -^s^^^,-^J.,5,5,. (1) 

The term — X](i 7) JijSiSj describes the random interaction between spins of the spin glass with gaussian-distributed 
random J^ of width J and the convention that («, j) counts pairs i,j only once and only with i ^ j. 
To investigate the physical effects which result from this model we express the partition function employing the replica- 
trick, perform the configuration average over the frustrated spin couplings, and after the decoupling procedure by 
use of the Hubbard-Stratonovich transformation obtain the saddle point equations for the model (for the detailed 
description of this technique see Ref .113) . The partition function of the model is thus obtained as 
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with /3 = 1/T. P{Jij) denotes the Gaussian distribution of couplings, the superscripts a relate to replicas. Here and 

in what follows the superscript z will be dropped for the hole-spin s. 

The Kondo-term — (k/A^i)s^jJj^ Si is decoupled with the help of the identity 
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The spin glass term in the exponent of the partition function is identical to the hamiltonian of the wellknown 
Sherrington-Kirkpatrick (SK) model (seella£3) . Thus the decoupling of this term is completely equivalent to that of 
the SK model. 
The whole decoupling procedure transforms the partition function into 
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where the fields Ra:VaTMa facilitate the decoupling of the Kondo-like interaction, while the field Q"^ does the job 
for the spin glass interaction. 
Now we define the replica-symmetric and spatially homogeneous saddle point Ansatz for the decoupling fields by 



Ra=V 

Ma ^ im 

Va = ip 

Q"^ = g, for a^h. 



(5) 
(6) 
(7) 
(8) 



For the fields Q°'^ , replica-diagonal and replica-offdiagonal components are distinguished in general. The replica- 
offdiagonal components are related to the SG order parameter, whereas the replica diagonal components are noncritical 
for the Ising spin glass modelMI. After the substitution of the saddle-point Ansatz in the partition function the free 
energy of the system is obtained as follows 
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The free energy becomes extremised by 77, p, m and q satisfying the following saddle point equations 

77 — p/Ni + m 



(9) 
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The value q — ^ {S°-S'')in.s + ^ j^^ ('5'"5'^)o»*« (a ^ h) is the weighted average of Edwards-Anderson order parameters 
inside and outside the polaron respectivelytS. p — (s") represents the thermodynamically averaged spin of the hole 



and m — l/A^i(X]j=i "^i") denotes the average magnetisation in the polaron. The equations ([T^)-([l4|) describe the 

whole phase diagram of the model in the static mean field approximation. From the qualitative analysis of these 

equations three characteristic temperature ranges can be identified. 

(i) At high temperatures T > k/\/Ni the equations have only trivial solutions. The system is paramagnetic and 

thermal fiuctuations suppress the magnetic polaron formation. 

(ii) At k/\/Ni > T > J the magnetic polaron is formed but the Mn system remains in the paramagnetic phase. 

(iii) At T < J the SG order is expected to show up in the Mn system. Hence the value n/ ^/Wl determines the 

temperature of the MP formation. In the experiments the MP is formed in the paramagnetic temperature regime, 

which implies the inequality k/\/Ni > J. 

III. TEMPERATURE DEPENDENCES OF MAGNETISATION IN THE POLARON AND OF THE SPIN 
GLASS ORDER PARAMETER IN THE VICINITY OF THE PHASE TRANSITION 

Here we examine the dependence of the SG order parameter q and of the cluster magnetisation m on T — J, where 
J corresponds to the transition temperature of an independent SG model. The main diffecence between the present 
case described by equation (nj) and the analogous equation for the independent spin glassta is caused by the presence 
of the exchange field 

P 

hexch ^ ^irr l^'-'j 

which the localised hole exerts on the pseudospin. It is worth noting that the influence of the exchange field on the 
magnetisation in the polaron and on the SG order parameter q is analogous to the infiuence of an external magnetic 
field localized in the volume of the magnetic polaron (see Eqs.(n3) and dl4)). The nonzero magnetic field competes 
with the exchange interaction between Mn spins and destroys partially their-ordering. This results in the smearing 
of the phase transition and lowering of the (smeared) transition temperaturell3. The transition should still show up 
in a change of the temperature behavior of the magnetization m and of the SG order parameter q. In order to see 
these expected changes we analyse the temperature behavior of the quantities m and q in the vicinity of the point 
T = J. The case when the polaron constitutes a fixed part of the Mn spin system is considered, i.e. a — ^ fixed, 



< a < 1. We expand the right-hand sides of Eqs.(|l3[), ( |14D for small q up to quadratic terms. After integration 
over the random field z the equation for q becomes 

2/3^*92 _^ (1 _ /32 j2(l - AC))q - C = 0, (16) 

where C = atanh {f3hexch)- The physical solution for q reads 



_ /3V2 - 1 - 4/3V^C + ^(1 - /3V2)2 ^ sp^J^C 
^~ 4/34 J4 • ^'^'> 

In the region |1 — P'^.P] ^ P'^J^C the behavior of q changes (for C = nonanalytically) at pj — 1. For T > T^ (i.e. 
in the paramagnetic phase) 

g«C«^^(l-20), (18) 
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where 9 = T/T^ - 1. In the SG-phase, T <Tc 
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q^2^^{l + 29)-e. (19) 

The dependence q{9) is linear both above and below T^, but one observes different slopes in the paramagnetic phase 
and in the SG phase. Similar behavior is observed in the magnetisation m. We find 

f %^ (1 - 9) for T > Tc, 

m^l u ^\^ , , (20) 

[ ^^ - a%f^ (1 + 2561) for T < Tc. ^ ' 

The linear dependence of m on becomes very weak below T^, i.e. in the SG phase. 

The solutions of the saddle point equations describe also the smearing of the otherwise sharp phase transition as 
caused by the finite exchange field in the magnetic polaron. The smearing is efficient around the point Tc = J within 
the temperature interval given by |1 — /3^J^| <C P^J^C. Beyond this interval the power law behavior of m and q 
characteristic of the transition in the pure spin glass is retrieved. 



IV. NUMERICAL SADDLE-POINT SOLUTIONS AND COMPARISON WITH EXPERIMENTS 



In order to obtain a detailed description of the physical behavior we solved the saddle-point equations (12)-(h4) 
numerically. It is convenient to define a composite SG order parameter according to g = aqi + (1 — a)q2, where 

qi= f dz tanh^ h{JV^z + K^) j , (21) 



(72 ^ / dz tanh-^ (PJy/qz) . (22) 

Jg 

describe spin glass order inside (qi) and outside (52) the magnetic polaron, respectively. 

A qualitative understanding of the coupled selfconsistent solutions involving polarisation p, magnetisation m, and 
the SG order parameters qi, ^2 can be obtained as follows. Fig.|l| displays the magnetisation, which is determined 
from Eq.(pl3[), and the SG order parameters. The thermal behavior of m and qi is strongly affected by the coupling 
Ki = k/Ni which controls the interaction energy of a pseudospin with the hole-spin. For a large value ki — S.lAmeV 
(curve 1 in FigJ^) the magnetisation m decreases monotonuously with temperature and the phase transition appears 
strongly smeared by a large exchange field (which keeps gi at high values due to non spin glass contributions and 
thus inhibits its critical decay at the zero field freezing temperature). At smaller ni — 0.628meV (curve 2 in Fig.n^) 
instead, m(T) displays a maximum near the almost sharp SG transition signalled by the rather steep decay of both 
SG order parameters gi, within the MP, and 92 • This is in accordance with the analytical solution given in Eq.(Effl). 
For fixed temperature, m and qi grow monotonuously with ki. The temperature dependence of the average hole-spin 
p at two different values of the interaction energy k between hole-spin and pseudospins within the MP is shown in 
Fig.0. The change of k affects the thermal behavior of the polarisation p{T) in a similarly strong way as a /«i-variation 
alters the magnetisation m{T). In turn, the spin glass order parameter 52 outside the magnetic polaron depends very 
weakly both on n and ni. Fig.n^ shows that with increasing ki the temperature dependence of gi, which resembles 
the g2~curve for small ni (almost sharp SG transition), changes over and approaches the magnetisation curve m{T) 
in the large Ki-regime. 

The total number N of pseudospins enters in all formulae only in the form A'^i — aN. Since, additionally, the 
distinction between magnetic order inside the MP and the one within the remaining spin glass domain (not exposed 
to the exchange field) will influence the quantitative agreement between theory and experiment, it is convenient to 
analyse thermal properties as a function of A'^i and a. Ni together with k determines the exchange energy ni on 
a pseudospin. Changing Ni, at fixed k and a, is thus equivalent to tuning ki, which was defined in the preceding 
paragraph. The parameter a determines the relative extension of the polaron and enters in the SG order parameter 
via the relation q = aqi + (1 — a) 52- Large values of a correspond to a strong influence of the polaronic exchange 
field on the spin glass order because a is the weight of the gi-contribution to the composite SG order parameter q. 
As a decreases, the peak m{T) shifts to higher temperatures and becomes more pronounced (see Fig.H). 
Measurements of tlie magnetic-field-induced circular polarisation degree of exciton luminescence (reflecting the mag- 
netic susceptibilityli3) and of the magnetic polaron energy (determined by selective excitation of localised excitonsEJ) 
against temperature are displayed in Fig.^. The polarisation degree shows a cusp at T = 8K characteristic of the 
spin glass phase transitionlid'EJ. The cusp temperature coincides with thCj-SG transition temperature, which was ob- 
tained for this manganese concentration by other measurements (see Ref.Elilj). The experiments on the temperature 
dependent EMP energy (solid circles in Fig. 3) do not show any peculiarity at the spin glass transition temperature 
Tc = 8K. However, a plateau is observed at lower temperatures. We argue that this plateau results from the onset of 
SG order inside the polaron. The plateau extends over a temperature range lying below the cusp temperature of the 
zero field spin glass, because of the large exchange field inside the MP. 

In order to determine the parameters of our model we fitted the experimentally known temperature dependence of the 
exciton magnetic polaron (EMP) energy in a Cdo.67Mno.33Te epilayer. The exchange part of this energy is calculated 
according to 

e = Kmp (23) 

We fix the value of J = 8K (the cusp temperature in the circular polarisation data) whereas the effective exchange 
interaction k between a hole spin and pseudospins (tetrahedrae of Mn spins) in the polaron is chosen as a fit parameter. 
We also vary the parameters N and a to obtain maximal agreement of the calculation with the experimental points. 
The best fit, given by the solid line in Fig.H, is obtained with a = 1, A^ = 20, k = AO.bmeV, which is already in 
a rather good agreement with the experimental data. This justifies our assumption that the magnetic contribution 



to the EMP energy dominates thermal variations and is hence responsible for the main features of its temperature 
dependence. The numerical calculations indicate that for these parameters the hole spin is completely polarised, 
which means p — 1 for all temperatures under consideration. Thus the temperature dependence of the EMP-energy 
reflects the polarisation degree of the Mn spins in the region of the excitonic wave function. Calculating the exchange 
field (~ K,/{aN)) in the polaron we obtain a vakic of about 4.7T, which is in reasonable agreement with the value 
5.9T obtained by other experimental techniqueal3. If we take into account that N is the number of psewrfospins 
representing tetrahedrae of Mn spins, we can as well estimate the number of Mn spins-inside the polaron to be about 
60, which is roughly a factor of 2 smaller than the value estimated from experimentQ'tj. 

Despite the rather good agreement between experiment and the model description discussed so far, there are two 
reasons why the present situation is not yet completely satisfactory. On one hand, a = 1 corresponds to a polaron 
extending through the whole spin glass, on the other, it is known that, besides the percolating spin glass cluster, 
there exist small clusters of antiferromagnetically coupled manganese spins. It turned out that these two problems 
are connected, since an extended model containing bulk spin glass order and clusters of a small number of antiferro- 
coupled Mn spins allows to describe the experimental data with a more realistic value of a. 

We model the temperature dependent polarisation of the small Mn-spin clusters by use of the wellkniawn modified 
Brillouin functiorO commonly used to fit the experiments on Zeeman splitting of excitons in (Cd,Mn)TeEl. This yields 

m., = xS^ffB,/, (1^) • (24) 

Here the effectpe spin Seff < 5/2 and the effective temperature T^ff reflect the antiferromagnetic coupling of Mn- 
spins (see Ref.EZI); x = 0.33 stands for the manganese concentration, and /i denotes the Bohr magneton. In the 
argument of the Brillouin function we insert instead of an external magnetic field the exchange field Bexch between 
hole spin and Mn-spins in the MP. For the manganese concentration of 33% the exchange field is Bexch ~ 5.9TE3. 
The magnetic part of the EMP energy is now given as the sum of the contribution originating in spin glass order and 
a second term representing the effect of small clusters; the total magnetic EMP energy then becomes 

Eemp = Kpm + Koprrici, (25) 

iVi + l\cl iVi + l\cl 

where Nd denotes the |Cumiber of spins within small clusters in the MP. The exchange constant kq — 880meV is 
known from experimentEJ. 

We fitted the experimental data with the form (P5|). Now, in addition to n, a and iVi, we face three fit parameters 
more. These are the number of spins in small clusters Nd and the parameters of the modified Brillouin function Seff 
and Teff . The aim of our procedure is now to prove that 

i) it is possible to fit experimental data with allowance for the additional presence of small clusters, keeping all 
parameter values within the physically meaningful regime, together with 

ii) the plateau in the temperature dependence as a genuine effect of the spin glass. This plateau cannot be obtained 
selfconsistently from the model of clusters of antiferromagnetically coupled spins only. The fitting procedure is 
illustrated by Fig.|[ 

As mentioned above the searched for function Eemp (T) is determined by the temperature dependences of the magneti- 
sations m and md because the hole-spin is completely polarised in the entire temperature region under consideration. 
For bulk-like real systems one expects a nonzero portion of the spin glass cluster to be out of reach of the exciton 
wavefunction, hence outside the polaron. Therefore the real physical situation corresponds to a < 1. The ratio 
Ni/Nd defines the relative contribution of the SG and of small clusters to the EMP energy. The lack of experimental 
information on the precise value of a and Ni/Nd forced us to analyse fits with different values of a for the interval 
from 0.3 up to 0.7 and for Ni/Nd of order 1. The results for the fits were remarkably insensitive to the variation of 
a at least within this regime. In the fit presented in Figs.H (dash-dotted line) and (solid line) we fix a = 0.5 and 
iVi — Nd, which we consider to be a good choice for the real system under consideration. While these relations were 
not rigorously dictated from experiment, the insensitivity of the fits with respect to variations of these two parameters 
over a wide range guarantees that the fit, given for the above choice, will describe the real system. 
Having chosen a and Ni/Nd, we are left effectively with four fit parameters, k, A^i, Seff , and Teff, which remain to be 
determined. In the following we define the four relations employed here to determine these parameters unambiguously. 
The first three of these conditions are the principal experimental features, which need to be reproduced by the fit 
curve: 

i) the absolute plateau-height at low temperatures (~ 28meV), 
ii) the high temperature end of the plateau ('plateau-ridge') at w iK, and 
iii) the decay of the magnetic EMP energy at higher temperatures, as seen in the experiment way above the spin glass 



freezing. 

The fourth condition is concluded from the foUowing consideration: The plateau ridge is strongly connected with the 

position of the peak caused by the SG contribution (dashed line in Fig.^), which itself is determined by the value 

Ki = k/Ni. Keeping the position of the peak as close as possible to the plateau ridge, condition iv) results and we 

obtain a single value ki for each k; hence the parameter iVi is now completely determined by the value of k. 

The matching of the plateau's high temperature end with the experimental data is achieved by tuning the form of 

the contribution of small clusters (by means of the parameters Tejf and Seff) for each k. This requirement provides 

a relation between k, Seff, and T^ff. 

The absolute value of the plateau is defined by the parameters Seff and k, whilst the role of Tg// and A^i is negligible. 

Thus the absolute plateau height of the energy curve yields a relation between k and Seff ■ 

Finally the form of the decay at high temperatures as seen in the experiments provides one more relation between 

the parameters k, Seff, Teff, and iVi. 

Fig.^ shows the form of the SG contribution (dashed line) and that of the small cluster contribution (dashed-dotted 

line) to the EMP energy together with the resulting combined curve of the temperature dependence (solid line) . Note 

that the form of the small cluster contribution is concave. It follows that the plateau results from the SG contribution 

to the EMP energy. 

One can try to fit the same experimental data with the modified Brillouin function alone. However, in this case one 

should take the parameters obtained from the Zeeman splitting fit in Cdo.67Mno.33Te, since now the Brillouin function 

describes the whole spin system and not only the part of small clusters as above. If one does so, the Brillouin fit docs 

not coincide with the experimental data. Agreement can only be achieved by using Brillouin function fit-parameters 

very different from those obtained by the Zeeman splitting. For this reason such a procedure does not result in a 

consistent description of the complex magnetic behavior of (Cd,Mn)Te. 

In conclusion, the fit described above shows that the main qualitative feature, which is the plateau in the temperature 

dependence of the EMP energy, is the specific result of the spin glass order in the polaron. The plateau occurs at a 

temperature lower than the bulk SG transition temperature because of the high exchange field in the polaron. 

V. DYNAMICAL ASPECTS OF MAGNETIC POLARON FORMATION 

In this section we suggest a phenomenological model for the dynamics of the MP formation. The approach presented 
here aims at taking into account the contribution of the SG-phase. Thus we do not consider the nonmagnetic 
contributions to the MP formation (prelocalisation) , the shrinking of the hole-wavefunction during MP~formation 
and [-t|lj LO || rn icroscopic mechanisms of the spin- and energy relaxation (a description of these subjects is given in 

We suppose the changes of the averaged spin of the hole and the magnetisation in the polaron to be slow in time 
in comparison with the thermodynamic fluctuations so that, at each instant of time, we can define the state of the 
system by thermodynamically averaged quantities (magnetisation, average hole-spin etc.). Then the formation of 
the MP can be considered as an evolution of MP parameters m and p in time from zero values at t = which is 
the trivial solution of the static mean field equations (|2| and |l^) towards the nontrivial static solutions at i ^ 00. 
This evolution is governed by an additional dynamical part. Hence we propose the following phenomenological set of 
equations 



i> + 7P + ri^p = Vl^ tanh [Pnm] 



(26) 
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This form of the dynamical part contains only the leading time derivatives and is linear in p and m. An alternative 

form containingonly first time derivatives does not yield a time evolution under the starting conditions p = 0, to = 0. 

The Eqs.(|2q), ( |27| ) describe the dynamics of nonlinearly coupled harmonic oscillators with damping. The characteristic 

frequency il is the oscillation frequency, while 7 characterises the damping. 

We assume q to satisfy equation ( p^ with time-dependent p and to,. This assumption is consistent with the definition 

of q as Edwards-Anderson order parameter. The equilibration time for q is of order l/T - the characteristic time 

of thermodynamic relaxation processes at this temperature. For the suggested model this time is shorter than the 

characteristic evolution time of p and m. 

The dynamics of the MP formation is described by the solution of equations (p^), (p7|), (111) with the initial conditions 

p{t = 0) = 0, m{t = 0) = 0. We can obtain an analytical solution for the initial stage of the MP formation where 



Pmn <C 1, Pnp/Ni <C 1 so that we may expand in these values up to the linear terms. In the paramagnetic phase q is 
nonzero only due to the exchange field of the MP. The relevant values of z in the integral of Eq. (gTJ) lie in the region 
J^/qz <^ Kp/Ni. Hence we neglect q in Eq.(|27|) and obtain for the initial step of MP-formation t <C min{l/j, 1/^) in 
the paramagnetic phase (T > J) 



p + JP + fl^p = n'^(3Km, 



(28) 



m + 7m + n m — il (3k 



Ni 



(29) 



The nontrivial solution of this system of linear differential equations with the initial conditions ?n(0) — 0,p(0) = is 

m = Ce"^*/^ sinh 



J-f^+m^ (Bk/^/nI 



1 



(30) 



Cv/iVie-T*/2 sinh 



J-i^+m^i^Pn/^i 



(31) 



One obtains a nonoscillating solution, which corresponds to the formation of the polaron only under the condition 
fin/ \/Ni — 1 > 0. The same inequality was obtained in the end of section ^ from the static selfconsistency equations 
as a necessary condition for the polaron formation. The constant C and the characteristic frequencies 7 and D, could 
be inferred from the comparison with experiment. 

Let us now consider the case of spin glass order. In the vicinity of the phase transition the integral in Eq. (E^ is 
dominated by /3Jy/qz ^ 1. Since /Snp/Ni <C 1, we have 



dz tanh 



PiJVQ^ + n^^) 



t'Anh{(3Kp/Ni)(l- / dztanh^iPJ^zU 
ta.nh{(3Kp/Ni){l -/(?)). 



(32) 



At times t <C I/7 the condition Jy/qz ^ np/Ni is fulfilled for the relevant domain of z in the integral. Then we can 
set I{q) « q with q evaluated from Eq.(|l4) with p = 0. Expanding tanh{x) « x the following equations are obtained 



p + jp + n p — ft (3Km, 



(33) 



m + jm + il m — il (3k — G{q), 
where G{q) = 1 — I{q) ^ I — q. The solutions of these equations read 

m = Ce^^*/"^ sinh 



^-^^^ + m^((3K^G{q)lN, 



(34) 



(35) 



p = C^/Ni/G{q)e-'''/^ sinh 



yj^ + An^(^pKVGiq)/N, 



(36) 



^From G ~ 1 — q < 1 it follows that the spin glass order slows down the magnetic polaron formation. 

The time dependence of the exchange part of the MP energy at the initial stage of MP formation (Eq.(^)) is given 

by 



£(i) ex V^i/G(g)e-^* sinh^ 



^j^+4n^(^(3K^G{q)/Ni~lj 



(37) 



where 



/-i( \ — j ^ i'^ paramagnetic phase ,„_^ 

^ ' ~ 1 1 — /(g) in spin-glass phase. ^ ^ 

Thus the MP formation depends exponentially on time. Representing the magnetic part of the MP energy as e(i) ex 
exp(i/rj), we obtain the following temperature dependence for the MP formation time 



Tf 



^72 + 4172 (/3«:v/G(g)/7Vi - l) - 7 j (39) 



The introduced model allows one to discriminate two limiting cases of the MP formation dynamics: the case of "strong 
dissipation" 7 ^ 51 and the case of "weak dissipation" 51 ^ 7. The expressions for the MP formation time in these 
two limit cases read 



Tf 



^(f3K^Giq)/N,-l)] \ for 7 » 51, 

^^ - (40) 

, for 51 > 7. 



2n {fiK^G{q)/Ni - 1 



Neglecting the contribution of the exchange field to q in the SG phase [q (x {J — T) / J) one can see the crossover 
of the temperature dependence of the MP formation time at the SG transition. At 7 ^ 5} the crossover is from 
T/ oc T in the paramagnetic phase to Tf ex T^/^ in the SG phase. At 51 ^ 7 the crossover is from r/ ex T^/^ in 
the paramagnetic phase to Tf c>c T^/"* in the SG-phase. The qualitatively different temperature behaviors of the MP 
formation time in the different extreme cases give a possibility to infer from the comparison with experiment the 
relevant corresponding d#imnical model. The limit case of "strong dissipation" coincides with the dynamical model 
introduced in the paperai3'&E3. 

Taking a typical value of the exciton magnetic polaron formation time of order lOOps at temperature T k, 21 
and the parameters k = ilmeV, Ni = 20 (the parameters of the solid line in FigJl), we can evaluate the order of 
magnitude of the characteristic frequences 51 and 7. For the "weak dissipation" case (51 ^ 7) we find 51 ^ lO^Hz. 
For the case of "strong dissipation" (7 ^ 51) we find 51^/7 ^ IQ^Hz, which implies 7 ^ 51 ^ IQ^Hz. The evaluation 
with the parameters k = 65m,eV and A^i — 50 corresponding to the dashed-dotted line in Fig.^ would give the same 
order of the characteristic frequencies. 

VI. CONCLUSION 

We elaborated a detailed model description of specific thermodynamic features of the magnetic polaron near the 
spin glass (SG) transition and due to SG order. The effect of the exchange interaction between carrier-spin and 
the spin glass is analogous to the effect of a spatially confined exchange magnetic field, which competes with SG 
order. Consequently, this field leads to both a decreased and smeared freezing temperature within the range of the 
magnetic polaron. The onset of spin glass order results in a qualitative change of the temperature dependence of the 
magnetic polaron energy as observed experimentally. The obtained solutions of the selfconsistent equations for the 
polaron magnetisation and for the thermodynamically averaged hole spin allowed to fit the experimental data of the 
temperature dependence of the exciton magnetic polaron energy in (Cd,Mn)Te. 

An extended model describing different dynamical scenarios of the magnetic polaron formation was introduced. These 
dynamical classes can be distinguished by specific thermal behavior of the MP formation time. A comparison with 
experiment on thermal behaviour should hence allow to identify the relevant dynamical scenario. 
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Figure 1 



FIG. 1. Temperature dependences of the magnetization m (solid lines) and SG order parameters qi (dashed hnes) and g2 
(dotted hnes) at two different values of exchange constant ki — k/N\: (a) ki = S.lAmeV and (b) ki — 0.628meV . Parameters 
K — 44meV and a — 0.7 are fixed. The curves are normalized to the largest value of the corresponding parameter. 
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Figure 2 



FIG. 2. The temperature dependence of the average spin of the hole p{T) for different k: (a) k = lOmeV, (b) k = 5meV, 
(c) K = 3meV. The parameters J = S-R', a = 0.7, h — 0, and ki = 0.63meV are fixed. 
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Figure 3 



FIG. 3. Temperature dependences of the magnetization m (solid lines) and the SG order parameter q — aqi + (1 — a)q2 
(dashed lines) at two different values of a: (a) a — 0.7 and (b) a = 0.2. Other parameters are fixed: J — 8K, ki = 0.628meV , 
K = 44meV . 
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Figure 4 



FIG. 4. The temperature dependence of BMP energy (closed circles) and the magnetic-field-induced polarization degree 
measured at 0.4T (open circles) of a MBE-grown epilayer CdQ.e7Mno.33Te. Results of theoretical calculations: Solid line: Fit 
with SG only, k = 41meV, N ^ 20, a = 1. Dashed-dot line: Fit with SG and smaU clusters, kq = 880meV, T^ff = 8K, 
Seff = 0.12, K = 65meV, N = 100, a = 0.5, A^i = N^i = 50. 
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Figure 5 



FIG. 5. The fitting procedure for the model of SG and small spin clusters. Solid circles: experimental points as in Fig.W. 
Dashed line: the SG contribution to the EMP energy. Dashed-dotted line: the contribution of small clusters to the EMP 
energy. Solid line: the resulting fit curve. The parameters of dashed, dashed-dotted and solid lines are the same as by the 
dashed-dotted line in Fig.^ 
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